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ABSTRACT 

We present an analysis of the observed broad iron line feature and putative 
warm absorber in the long 2001 XMM-Newton observation of the Seyfert-1.2 
galaxy MCG-6-30-15. The new kerrdisk model we have designed for simulat- 
ing line emission from accretion disk systems allows black hole spin to be a free 
parameter in the fit, enabling the user to formally constrain the angular mo- 
mentum of a black hole, among other physical parameters of the system. In an 
important extension of previous work, we derive constraints on the black hole 
spin in MCG-6-30-15 using a self-consistent model for X-ray reflection from the 
surface of the accretion disk while simultaneously accounting for absorption by 
dusty photoionized material along the line of sight (the warm absorber). Even 
including these complications, the XMM-Newton/EPlC-pn data require extreme 
relativistic broadening of the X-ray reflection spectrum; assuming no emission 
from within the radius of marginal stability, we derive a formal constraint on the 
dimensionless black hole spin parameter of a = 0.989io;QQ2 a ^ 90% confidence. 
The principal unmodeled effect that can significantly reduce the inferred black 
hole spin is powerful emission from within the radius of marginal stability. Al- 
though significant theoretical developments are required to fully understand this 
region, we argue that the need for a rapidly spinning black hole is robust to phys- 
ically plausible levels of emission from within the radius of marginal stability. In 
particular, we show that a non-rotating black hole is strongly ruled out. 

Subject headings: accretion, accretion disks - black holes - galaxies:nuclei - X- 
rays:spectra 



1. Introduction 

Accreting black holes are the driving force behind some of the most powerful processes 
in the universe. The central regions of active galactic nuclei ( AGN) , in particular, are prodi- 
gious sources across the electromagnetic spectrum including the X-ray band. The X-ray 
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continuum (which is accurately approximated by a power-law up to ~ 100 keV) is thought 
to be produced by Comptonization processes in the corona (or base of a jet) surrounding 
the inner part of the accretion disk. A portion of the X-ray photons produced in the corona 
irradiate the underlying optically-thick disk, producing the so-called "X-ray reflection" sig- 
natures in the observed spectrum (Guilbert k Rees 1988; Lightman k White 1988). These 
X-ray reflection signatures consist of fluorescent and recombination emission lines sitting on 
a continuum due to Compton scattering and the summed radiative recombination continua 
of the excited ions in the photoionized disk surface (George k Fabian 1991; Ross k Fabian 
2005). The Fe-Ka line is the most prominent of these features due to its energy (at 6.4 keV 
it is visible above the direct power law continuum), and the high astrophysical abundance 
and fluorescent yield of iron. Additionally, this line is often significantly broadened by both 
the standard Doppler effect and relativistic processes, the effects of which increase the closer 
the line is emitted to the event horizon (Fabian et al. 1989; Laor 1991). These relativistic 
processes include light bending, beaming and gravitational redshifting, all of which result 
in a greatly elongated and skewed line profile; in particular, the line profile can display an 
extended low-energy tail primarily resulting from gravitational redshift. The line emission 
region is sufficiently close to the black hole that frame-dragging effects associated with the 
black hole spin can be important in determining the line profiles. The broad iron line is there- 
fore a powerful probe of the relativistic effects on the spacetime immediately surrounding 
the black hole. 

The first broad iron line robustly detected and resolved in an AGN was found by the 
Advanced Satellite for Cosmology and Astrophysics (ASCA) in the Seyfert-1.2 galaxy MCG- 
6-30-15 (Tanaka et al. 1995; Iwasawa et al. 1996), and since then has been extensively stud- 
ied with BeppoSAX (Guainazzi et al. 1999), RXTE (Lee et al. 1999, 2000), Chandra (Lee 
et al.2002; Young et al. 2005) and XMM-Newton (Wilms et al.2001; Fabian et al. 2002; 
Reynolds et al. 2004). All of these results show that the broad iron line feature is consistent 
with a highly redshifted line from the inner parts of an accretion disk; no alternative hy- 
pothesis has yet explained the spectrum of MCG-6-30-15 satisfactorily (Fabian et al. 1995; 
Reynolds k Wilms 2000; Vaughan k Fabian 2004; Young et al.2005). Subsequent ASCA, 
Chandra and XMM-Newton studies have discovered similarly broad iron line profiles in sev- 
eral other Seyferts, such as MCG-5-23-16 (Dewangen, Griffiths k Schurch 2003), NGC 3516 
(Turner et al.2002), Mrk335 (Gondoin et al.2002), and Mrk766 (Pounds et al. 2003a). 

In addition to having a sample of objects that have been observed with robust broad 
iron lines, it is equally important to have a precise model to use in fitting the data. The 
two line profiles currently included as standard in the current versions of the commonly 
used spectral fitting package xspec (Arnaud 1996) are useful as a starting point, but ul- 
timately quite limited in terms of their ability to accurately parameterize the line. The 
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diskline model (Fabian et al. 1989) describes the line profile from a disk around a non- 
rotating (Schwarzschild) black hole, and, due to the approximations employed, does not 
include relativistic light bending. Similarly, the laor model (Laor 1991) has important con- 
straints as well: this is a fully-relativistic model, but the dimensionless spin parameter of 
the black hole is hard-wired at a = 0.998, the equilibrium spin of a black hole accreting 
from a Novikov & Thorne (1974) accretion disk. Furthermore, due to the computational 
realities of the early 1990s, the relativistic transfer functions underlying the laor model are 
pre-calculated and tabulated, yielding noise (or even gross inaccuracies) in the line profiles 
produced, especially at very high disk inclination angles. Given these limitations, as well as 
the high quality of AGN spectra currently being obtained with Chandra and XMM-Newton 
(and hopefully in the future with the remaining operable instruments aboard Suzaku), it is 
imperative that X-ray astronomers have access to effective models that are fully relativistic, 
accurate, and that allow black hole spin to be fit as a free parameter. This latter point is 
crucial if we are to hope to constrain the spin of astrophysical black holes using broad iron 
lines and other X-ray reflection signatures. 

Three new relativistic line models have recently been developed for this purpose and im- 
plemented in a form that can be readily used by X-ray astronomers: the ky suite (Dovciak 
et al. 2004) and similar codes by Beckwith & Done (2004) and Cadez & Calvani (2005). 
These models achieve comparable results for the morphologies of the line profiles, and all 
offer significant improvements over the diskline and laor results in terms of accuracy and 
precision over a wider range of physical parameters. Most importantly, these models leave 
the spin of the black hole as a free parameter, include emission from within the innermost 
stable circular orbit in the disk, and compute fully relativistic photon transfer functions. 
In practice, the xspec modules that implement both the ky and Beckwith & Done models 
use relativistic transfer functions stored in very large (multi-gigabyte) pre-calculated tables. 
Using these models, both groups conclude that fitting broad iron lines cannot truly constrain 
the spin of the black hole. Since their models allow emission from any radius of the accre- 
tion disk outside of the event horizon, they can produce iron line profiles with arbitrarily 
redshifted wings even if the underlying black hole spacetime has no spin. One must consider 
the physicality of this assumption, however: it is not possible to get line emission from radii 
deep within the plunging region for several reasons. Firstly, the closer one gets to the event 
horizon, the smaller the geometrical area of the disk available for emission becomes, while 
the gravitational redshift affecting any emission from such regions increases correspondingly. 
Both of these effects would tend to minimize the contribution to the line profile from most of 
the emission inside the radius of marginal stability. Further, the ionization of the disk ma- 
terial within ~ 4 r g of a Schwarzchild black hole is inevitably too high to produce significant 
line emission in the first place, and the optical depth of this material also decreases rapidly 
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within the radius of marginal stability (Reynolds & Begelman 1997; Young et al. 1998). For 
all of these reasons, one simply cannot produce significant line emission from any arbitrary 
radius outside the event horizon. This must be taken into account when modeling such line 
emission. 

In this paper, we undertake a rigorous re-analysis of the broad iron line seen in the 
~ 350 ks July-2001 XMM-Newton observation (Fabian et al. 2002). We assess the constraints 
on the iron line profile and black hole spin including the complications introduced by other 
spectral components displayed by this system, especially the substantial column of absorbing 
photoionized gas seen along the line of sight to the central disk of this AGN. To facilitate 
this investigation, we present a new variable-spin accretion disk line profile model, kerrdisk, 
that we are currently developing for public use in the xspec package. The principal dif- 
ference between kerrdisk and the previously mentioned variable-spin line profile models is 
that our model does not rely on gigabyte-sized, precomputed tables of the photon transfer 
function over many thousands of radii within the disk. Rather, as we will discuss below, 
we use the approach of Cunningham (1975) to decompose the line profile integral into a 
simple analytic component (describing the effects of relativistic beaming) and an additional 
transfer function (describing the effects of light bending) which is a slowly varying function 
of black hole spin and disk inclination angle and hence only needs to be coarsely tabulated. 
This technique allows us to obtain high-quality line profiles with very modest sized (40 MB) 
tables, making it substantially more portable and convenient than some other models. This 
accuracy, portability and comparative speed of computation make kerrdisk a valuable ad- 
dition to the X-ray astronomer's software arsenal. Additionally, a trivial change of the driver 
code allows the user to implement any desired limb-darkening law or compute the line pro- 
files to arbitrary resolution and accuracy (if, for example, the user is interested in the narrow 
cusps of the line profile). Neither of these modifications is possible if the limb-darkening and 
resolution are hard-wired into the tabulated relativistic transfer functions. 

In §2, we motivate our study by giving a brief recap of the history of the broad iron 
line in MCG-6-30-15. §3 describes our new iron line model, including a detailed comparison 
of kerrdisk to the other line profile models already discussed. §4 presents our analysis 
of MCG-6-30-15 using these variable-spin line models. This work is performed under the 
"standard" assumption that there is no observable X-ray reflection from inside the radius 
of marginal stability of the accretion disk. A future paper will extend the analysis to the 
more general case where physically reasonable levels of emission from this region are allowed. 
Under this assumption we arrive at the conclusion that the black hole in MCG-6-30-15 is a 
rapid-rotator. §5 and §6 discuss our results and draw conclusions, respectively. 
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2. A Brief History of Broad Iron Line Studies in MCG-6-30-15 

MCG-6-30-15 is an SO-type galaxy in the constellation of Centaurus that hosts a Seyfert- 
1.2 nucleus. It has a measured redshift of z = 0.008, placing it at a distance of d ~ 37Mpc 
using WMAP cosmological parameters. X-ray studies of this AGN have revealed a powerful 
central source along with a significantly broadened Fe-Ka feature. This broad line has 
led MCG-6-30-15 to become one of the most studied AGN in the X-ray band due to its 
potential as a probe of black hole and accretion disk physics. Although the mass of the hole 
is not yet well constrained, estimates based on X-ray studies have placed it in the range 
of 10 6 - 2 x 1O 7 M (Nowak & Chiang 2000; Reynolds 2000). McHardy et al. (2005) have 
further narrowed this mass range to 3 — 6 x 10 6 M Q using a variety of methods such as the 
M-cr relation, line widths from optical spectra, BLR photoionization arguments and X-ray 
variability. 

Reflection signatures from neutral material had been observed with EXOSAT (Nandra 
et al. 1989) and Ginga (Nandra, Pounds & Stewart 1990; Matsuoka et al. 1990), but it was 
not until deep observations of MCG-6-30-15 by ASCA that the broadened and skewed iron 
line was robustly detected (Tanaka et al. 1995). It was determined that the line profile 
matched that expected due to X-ray reflection from the surface of a relativistic accretion 
disk; the robustness of this interpretation was demonstrated in Fabian et al. (1995). Broad- 
band BeppoSAX data taken by Guainazzi et al. (1999) confirmed the Tanaka et al. detection 
of a broadened and redshifted iron line with an equivalent width of EW 200 eV. 

A detailed re-analysis of this observation by Iwasawa et al. (1996) identified a period 
of time when the source entered the so-called "deep minimum" state, marked by low con- 
tinuum emission. While in this state, the iron line width (and especially the extent of the 
red-wing of the line) markedly increased to the point that emission from a disk around a 
Schwarzschild black hole truncated at the radius of marginal stability (r ms = 6 r g ) could no 
longer reproduce the observed line profile. Noting that the radius of marginal stability moves 
inwards (to a location characterized by higher gravitational redshift) as the black hole spin is 
increased, it was subsequently argued that the black hole in MCG-6-30-15 must be rapidly 
rotating. Fitting sequences of Novikov & Thorne (1974) models to the deep minimum ASCA 
data (and assuming that the X-ray irradiation tracks the local dissipation in the underlying 
disk), Dabrowski et al. (1997) derived a lower limit of a > 0.94 on the rotation of the black 
hole. At this point, however, Reynolds & Begelman (1997) noted that physically plausible 
scenarios could result in sufficient X-ray reprocessing (including ionized iron line emission) 
from within the radius of marginal stability to fit the ASCA deep minimum state with even 
a Schwarzschild black hole. In order to illuminate the region of the disk within r ms , they hy- 
pothesized an X-ray source on the symmetry axis some height above the disk, and suggested 
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that iron line profile changes (and some part of the X-ray continuum flux changes) could be 
attributed simply to changes in the height of this X-ray source. 

After the discovery of its broad iron line feature, MCG-6-30-15 became the subject of 
many more observations. With the first observations of XMM-Newton in 2000, astronomers 
had a new tool with unparalleled throughput with which to examine this source in finer detail. 
The first XMM-Newton observation of MCG-6-30-15 was fortunate enough to catch the 
source in its "deep minimum" state characterized by low continuum flux and a broader than 
normal Fe-Ka profile. The resulting high signal-to-noise spectrum revealed an extremely 
extended red-wing to the line profile extending down to 3—4 keV (Wilms et al. 2001). Because 
so much emission seemed to be coming from radii deep within the gravitational potential well, 
and because the emissivity index was correspondingly quite high, the authors hypothesized 
an interaction between the spinning black hole and its accretion disk. Magnetic torquing 
effects between the two could result in the extraction of rotational energy from the black 
hole which would, in turn, power the high coronal emission seen from the inner radii of the 
disk (Agol & Krolik 2000; Garofalo & Reynolds 2005). An alternative mechanism consists 
of the strong gravitational focusing of a high-latitude source above a rapidly rotating black 
hole (Martocchia & Matt 1996; Miniutti & Fabian 2004). By contrast with the ASCA era, a 
Schwarzschild model appeared not to work for these data — a line extending down to ~ 3 keV 
would require line emission from extremely deep within the plunge region: r pa 3 r g = 0.5 r ms , 
a situation which appears unphysical due to the high ionization expected in this part of the 
flow (Reynolds & Begelman 1997). 

Reynolds et al. (2004) performed a follow-up and more detailed analysis of the XMM- 
Newton observation of MCG-6-30-15 taken by Wilms et al. (2001). They explicitly demon- 
strated that the iron line profile was inconsistent with an X-ray irradiation profile that follows 
a Novikov & Thorne (1974) dissipation law, even for an extremal Kerr black hole. Employing 
the generalized thin disk model of Agol & Krolik (2000) that includes a torque applied at 
f = r ms , Reynolds et al. (2004) suggested that this torque has a great impact on the disk 
emission seen in this observation. In the "deep minimum" state, the torqued-disk scenario 
requires that the disk is largely emitting via the extraction of rotational energy from the 
black hole rather than via accretion. These authors also examine spectral variability dur- 
ing the deep minimum state. By examining difference spectra and direct spectral fits to 
10 ks segments of data from this observation, the authors concluded that the intensity of the 
broad line seems to be proportional to the hard 2 — lOkeV flux of the source: the equivalent 
width of the line remains approximately constant while the source fluctuates substantially 
in amplitude. Such behavior is consistent with simple X-ray reflection models. 

The longest XMM-Newton observations of the source to date was published by Fabian 
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et al. (2002). This group found MCG-6-30-15 in its normal state, and recorded data for over 
87 hours. In this state, the bulk of the iron line emission was in a narrower line compared 
with the "deep minimum" state, although a very extended red-wing was still evident. The 
time-averaged EPIC-pn spectrum again showed that the Fe-Ka feature was very strong, and 
the long data set was of sufficient resolution and quality that the spectrum demanded a fit 
incorporating a full reflection model. Relativistic smearing needed to be applied not just to 
the cold iron line, but to the entire reflection continuum. Taking this into account, the iron 
line was once again found to produce emission within the radius of marginal stability for a 
Schwarzchild black hole. Fitting the line with the maximally spinning (a = 0.998) black hole 
laor model suggested that the line emissivity followed a e oc r~ a dependence with a = 4.5 — 6 
between an inner radius r m i n < 2r g and a break radius « 6r g . Beyond the break radius, 
the emissivity profile flattened to a ~ 2.5. This broken power-law form was strongly preferred 
by the data over the usual simple power-law emissivity functions usually fitted to such data 
- again, this reflects the high quality of the data. Fabian et al. (2002) also note that, in its 
normal state, difference spectra of MCG-6-30-15 show spectral variability from 2 — lOkeV 
in the form of a power-law: the iron line flux changed little between successive 10 ks frames 
of the observation, whereas the continuum flux varied by as much as a factor of ~ 2. This 
result is in contrast to that found by Reynolds et al. (2004), who observed an iron line flux 
proportional to the 2 — 10 keV continuum flux when MCG-6-30-15 was in its "deep minimum" 
state. Comparing these two studies, it appears that the iron line flux is proportional to the 
observed X-ray power-law continuum at low fluxes and then "saturates" to an approximately 
constant level once the observed X-ray continuum exceeds a certain level. This complexity 
could be due to light-bending effects if the power-law X-ray continuum source is situated 
close to the spin axis of the black hole (Miniutti & Fabian 2004). Alternatively, patchy 
ionization of the disk surface might produce such a saturation (Reynolds 2000). 

Dovciak et al. (2004) utilize their ky suite of iron line profile models to fit the time- 
averaged spectrum of MCG-6-30-15 as well, also using data from the long XMM-Newton 
observation of Fabian et al. (2002). They fit the 3 — 10 keV spectrum with four combinations 
of ky models, all involving a broad Fe-Ka kyrline, a narrow Gaussian emission line at 
6.9 keV (likely an ionized line of iron, as cited in Fabian et al. (2002)), and a Compton- 
reflection continuum from a relativistic disk (smeared using a kyconv kernel). The authors 
also find that models describing the disk emissivity as a broken power-law rather than a 
single power-law in radius achieve significantly better statistical fits, though among these 
broken power-law models, comparable fits can be obtained for a wide variety of black hole 
spin values. These authors therefore conclude that iron line profiles are not a good way 
to constrain black hole spin. At some level, this objection amounts to obvious point that 
arbitrarily large redshifts can be obtained around any black hole (rotating or not) if one is at 
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liberty to produce emission from any radius arbitrarily close to the event horizon. One must, 
however, consider the physical realism of this assumption in general, and the best-fit spectral 
parameters determined by these authors in particular. For their low-spin (a = 0.25) best- 
fitting model, Dovciak et al. concluded that a substantial amount of iron line emission had 
to originate from deep within the radius of marginal stability; translating the fit parameters 
of Dovciak et al. into standard quantities, they require an inner emissivity profile of e oc r~ 9 
starting at an inner radius of r in = 3.2 GM/c 2 , to be compared with the radius of marginal 
stability for a a = 0.25 black hole which is at r ms « 5.2 GM/c 2 . As has already been noted 
and will be explored further in §4, it is hard to understand how this region of the accretion 
flow could contribute to any part of the observed iron emission given that the fact it will be 
very tenuous and extremely highly photoionized (Reynolds & Begelman 1997; Young, Ross 
& Fabian 1998). 

The medium resolution (CCD) data discussed so far leave ambiguous the possible role of 
complex ionized absorption in distorting the observed X-ray continuum shape and hence the 
inferred iron line profile. To assess the role of this absorption, Young et al. (2005) performed 
and analyzed a deep, 522 ks grating observation of MCG-6-30-15 in May-2004 with the 
HETGS instrument on Chandra. This observation produced two important results: first, 
the authors found that the difference in the hard continuum spectrum between the high and 
low flux states was well described by a power-law of photon index T = 2.0^q 2 . This finding 
agreed with previous studies that indicated that the spectral variability of MCG-6-30-15 in 
its normal state is dominated by a power-law component (e.g. Fabian et al. 2002). Second, 
and most importantly, ionized absorption models whose continuum curvature mimics the 
red-wing of a broad iron line from 3 — 6keV were ruled out. Such models generically predict 
strong Ka absorption lines of intermediately ionized iron. Young et al. showed that these lines 
are conclusively absent, falsifying the ionized absorption model and further strengthening 
the relativistic smearing hypothesis. 



3. The kerrdisk Model 

3.1. Model Description 

As discussed in the introduction, we are motivated to construct an iron line model that 
allows the spin parameter of the black hole to be a free parameter and is maximally flexible. 

Our new relativistic emission line code, kerrdisk, is written in for.tr.an77 so that it 
can be easily meshed with xspec, and has now been compiled successfully on both Solaris 
and Linux platforms. The dimensionless spin parameter of the black hole (a = cJ/GM 2 ; 
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where J is the angular momentum of a black hole of mass M) can take on any value in the 
range — 1 < a < 1, where negative values of a correspond to a black hole that is rotating 
in a retrograde sense relative to the accretion disk. Although, according to the equations of 
General Relativity, a could have any arbitrary value, the cosmic censorship hypothesis states 
that naked singularities cannot exist in the universe, so a black hole must be shrouded by 
an event horizon. This limits the acceptable range of spin parameters to — 1 < a < 1. For 
simplicity we consider only prograde spins up to the Thorne (1974) spin-equilibrium limit, 
i.e., < a < 0.998. 

The limiting value of a = 0.998 for black hole spins was first discussed by Thorne (1974). 
Therein, the author contends that if black holes were simply accreting matter their spins 
would ascend to a ~ 1 rather quickly, but because material in the accretion disk radiates, and 
some of that emitted radiation is swallowed by the hole, a counteracting torque is produced. 
The origin of this counteracting torque lies in the photon capture cross-section of the hole: 
black holes have a higher capture cross-section for photons of negative angular momentum 
(opposite to that of the hole itself) than for photons of positive angular momentum. Thus, 
the accretion of emitted radiation from the disk results in an overall reduction in the angular 
momentum of the hole until it reaches a theoretical equilibrium value of a = 0.998. Recent 
work on magnetohydrodynamic (MHD) accretion disks suggest that the continued transport 
of angular momentum from matter within the radius of marginal stability, as well as angular 
momentum lost from the rotating black hole itself via Blandford-Zjanek like mechanisms 
(Blandford & Znajek 1977), may lead to a rather lower equilibrium spin (e.g., see the general 
relativistic MHD simulations of Krolik, Hawley & Hirose (2005)). Equilibrium spins as low 
as a ~ 0.90 are within the realm of possibility. 

Following the method of Cunningham (1975), we compute the line profile by employing 
a relativistic transfer function. For a given black hole spin a, disk inclination i, the observed 
accretion disk spectrum can be written as 

F ohs (E ) oc / =I(E o /g,0 e ,r e )ifj(r e ,g*;a,i)dg*r c dr e . (1) 

J yjg*(l-g*) 

Here, I(E,6 e ,r e ) is the rest-frame specific intensity of the disk at radius r e and energy E 
emitted at a angle 9 C to the disk normal, ip is the Cunningham transfer function which 
accounts for the effects of light-bending on the observed solid angle from each part of the 
disk, g is the ratio of the observed to the emitted photon energy (g = E Q /E em ), and g* is 
the relative value of g with respect to the redshift extremes obtained from that annulus, 

* _ 9 ~ ^min^c) /„n 

( \ — ( V ^ ' 

This formulation is attractive from the computational point of view; light bending effects 
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are isolated from beaming effects and encoded in the transfer function ip(r e , g*; a,i) which 
is a slowly varying function of its variables. This allows us to perform full computations of 
the transfer function at a relatively sparse set of points in parameter space, and then use 
interpolation to accurately determine the transfer function at a general point. 

Specializing to the case of a 5-function emission line allows us to eliminate the g* inte- 
gration from Eqn. (1). Assuming that the surface emissivity of the emission line as a function 
of emission radius r e and angle 9 e is f(r e ,9 e ), this gives 



We use the algorithms of Speith et al. (1995) to compute the transfer function i/j(r e , g*; a, i) 
as well as g m i n (r e ) and g m ax( r e), thereby allowing us to perform this integration. For a given 
spin a and inclination i, high quality line profiles (and all of the line profiles presented in this 
paper) are produced using a transfer function computed at 50 radial bins, equally spaced in 
the variable and 20 linearly spaced relative redshift (g*) bins across the line profile 

at any given radius. 

These values were determined via rigorous trial and error to make the best, smoothest 
possible line profile while keeping computation time to a reasonable length. Increasing either 
the number of radial or relative redshift bins did not significantly improve the integrity of 
the line profiles produced. Also, because the transfer function varies quite slowly with radius 
and relative redshift, greater frequency of sampling was not warranted. 

Our choice for the radial spacing corresponds to equal spacing of non-relativistic Kep- 
lerian velocities. We then evaluate the line profile integral (Eqn. 3) with a densely spaced 
logarithmic grid in r e , using linear interpolation in r e of the (rather sparsely sampled but 
slowly varying) transfer function. In practice, we use 5n e -zones per decade of r e where n e is 
the number of frequency elements in the output line profile (specified by the response matrix 
of the data being fit). Experimentation shows that this produces high quality line profiles 
even at very high resolution (i.e., large n e ). 

The Cunningham transfer function ip(r e , g*; a, i) is computed for a 20 x 20 grid in (a, €)- 
space and is stored in a 40 MB table that is accessed by the driver script for kerrdisk. 
Again, due to the slowly varying nature of the transfer function, simple linear interpolation 
to arbitrary spins and inclination angles produces accurate and high quality line profiles. 

For all of the models explored in this paper, we follow Fabian et al. (2002) and assume 
a line emissivity characterized by a broken power-law between some inner radius r min and 



F obs (E ) oc 




yjg*{l — g*) flWax — <?min 



ip(r e ,g*;a,i) f(r e , 9 e ) r e dr e . 



(3) 
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outer radius r max , i.e., 







fire) = < 



(r c /r br )- ai 
{r c /r hr )- a2 







'"e *^ ^min 
'"min ^ '"e *^ '"br 
'"br ^ '"e *^ '"ma> 
'"e — ^max 



(4) 



We note that, at present, the publicly available Speith algorithms on which this work is 
based do not support the proper computation of emission from within the radius of marginal 
stability, i.e., we are restricted to r min > r ms . Hence, all kerrdisk line profiles currently 
assume that the emission profile is truncated within r = r ms . 

The code takes ~ 4 s to produce a single line profile on a 2 GHz processor linux machine. 
When the model is used to fit data within xspec, several hundreds or thousands of iterations 
of the code are often needed to get a good fit, resulting in run times on the order of an hour or 
more. In principle, this run time can be significantly reduced by employing a parallelization 
scheme when fitting. Such a parallelization scheme is possible to implement in the ISIS 
spectral analysis package (Houck 2002) using the PVM module (Michael Nowak and Andy 
Young, private communication). 

The best way to examine the results of the code line profile simulations is to look at how 
the morphology of the line changes with alterations made in certain important variables. For 
illustrative purposes, Fig. 1 shows model line profiles with a rest-frame energy of E = 1 keV. 
The line profiles shown in this figure have the same parameter settings as those shown 
in the review by Reynolds & Nowak (2003) to enable easy comparison. The Reynolds & 
Nowak profiles are also based on the Speith algorithms but sample the transfer function 
more sparsely and employ a cruder integration technique to evaluate the line profile. The 
fact that our line profiles agree with, but are much less noisy than the Reynolds & Nowak 
profiles validates our newer integration/interpolation technique. 
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Fig. I. — Variation of the kerrdisk line profile with (a) disk inclination angle, (b) disk 
emissivity index, and (c) black hole spin parameter. In each case the range of emission in 
the disk is from r = r ms — 50 r g . In (a), a = 0.998, a = 1.5. In (b), a = 0.5, % = 40. In (c), 
i = 40, a = 3. 



3.2. Comparison With Other Relativistic Disk Line Models 



Verification of our new line profile code can be demonstrated through detailed compar- 
isons with existing public models, such as laor and diskline (where a = 0.998 and 0.0, 
respectively), as well as those of Dovciak et al. (2004) and Beckwith & Done (2004). Since 
the kyrline model of Dovciak has already been compared with the kdline model of Beck- 
with & Done (2005), we need only to compare kerrdisk with kyrline. These comparisons 
are shown in Figs. 2-5 and discussed in this Section. 
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Fig. 2. — Various iron line models for a Schwarzchild black hole. The inclination angles 
represented are 5° for (a), 45° for (b), and 80° for (c). Here a± = a% = 3.0, and r min and 
r mffl are held constant at 6r g and 50 r g , respectively. The diskline profile is in solid black, 
the kerrdisk profile is in dashed red, the kyrline profile including emission from within 
the ISCO is in dash-dotted green and the kyrline profile not including ISCO radiation is 
in dotted blue. 
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Fig. 3. — Various iron line models for a maximally spinning Kerr black hole. The inclination 
angles vary as above for the Schwarzchild case in (a)-(c). Other parameters are the same as 
those used in the Schwarzchild case, but now r min = 1.235 r g , corresponding to the radius 
of marginal stability for a maximal Kerr black hole, rather than the 6r g Schwarzchild r ms . 
The laor profile is in solid black and the rest of the color scheme is the same as that used 
in Fig. 2 above. 
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Given the lack of relativistic light bending in the diskline model, we expect our (fully 
relativistic) kerrdisk line profiles with a = to differ slightly from those computed with 
diskline, especially at large inclination angles. We indeed see slight differences (Fig. 2). 
We should stress, however, that the differences are minor and diskline should still be 
considered a perfectly acceptable model for disks around Schwarzschild black holes for all 
but the highest signal-to-noise data. Examining the comparison of the laor models with 
the a = 0.998 kerrdisk model, we note good agreement (Fig. 3). The slight differences 
that do exist are caused by an artificial smoothing of the laor line due to interpolation 
of a sparsely sampled transfer function. Again, however, we note that the laor model is 
perfectly acceptable model for disks around a = 0.998 black holes for all but the highest 
signal-to-noise data. 

The real power of this new generation of line profile models is the freedom in the spin 
parameter of the black hole, therefore the real verification of our code lies in a detailed 
comparison of kerrdisk and kyrline. Indeed, as shown in Figs. 2-3, the line profiles for 
kerrdisk and kyrline are virtually indistinguishable when one does not include radiation 
from within the radius of marginal stability in kyrline. Also shown in these figures are 
the effects of including emission from within the radius of marginal stability (down to the 
horizon, for the sake of illustration), as computed by Dovciak et al. (2004). This has the 
greatest effect on the line profiles from the slowly spinning holes since it is only in these 
systems that an appreciable fraction of the plunging region is subject to modest (rather 
than extremely large) gravitational redshift. We intend to extend our model to include the 
region within the radius of marginal stability in future work, carefully considering the effects 
of both the level of ionization and the optical depth of the material within the plunging 
region, in particular, as discussed in §1. Both of these effects can have a substantial impact 
on the contribution of emission from this region to the overall line profile. 



3.3. The Convolution Model 

In reality, the irradiation of the disk by the primary X-ray source results in a whole 
X-ray reflection spectrum consisting of Compton and radiative recombination continua plus 
numerous fluorescent and radiative recombination lines. The iron line is the most prominent 
due to its high rest-frame energy and its intrinsic strength, but many other species are also 
excited (e.g. oxygen, nitrogen, silicon and sulfur). The whole X-ray reflection spectrum from 
the disk will be subject to the same extreme Doppler and relativistic processes that combine 
to alter the morphology of the iron line. The most physically realistic simulation of the 
reflected spectrum from a photoionized disk surface has been presented by Ross & Fabian 
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(2005) - - these high-quality models capture the "traditional" X-ray reflection processes 
(Compton scattering, photoelectric absorption and fluorescent line emission) as well as the 
powerful soft X-ray radiative-recombination line emission expected from an X-ray irradiated 
photoionized surface of an optically-thick accretion disk. Ross & Fabian (2005) have provided 
their results in the form of tabulated spectra that can be used in xspec. We then wish to 
convolve this spectrum in velocity space with a relativistic smearing kernel such as the 
one used to generate the kerrdisk line profile. To facilitate this, we have also produced 
a convolution form of our line profile model, kerrconv, whose results mirror those of the 
kyconv model designed by Dovciak et al. (2004) and the kdconv model of Beckwith & Done 
(2004), provided that no radiation from within the ISCO is included in the latter two models, 
as was the case for the comparison of kerrdisk with kyrline and kdline in §3.2 above. 
As with kerrdisk, the kerrconv parameters include emissivity indices for the inner and 
outer disk separated by a break radius, inner and outer radii for the disk emission, the spin 
parameter of the black hole, and the inclination angle of the disk with respect to our line of 
sight. The kerrconv model is readily validated by applying it to a narrow emission line and 
comparing the resulting spectrum with the regular kerrdisk model. 



4. Determining the Spin of the Black Hole in MCG-6-30-15 

In this Section, we use our new models (kerrdisk and kerrconv) along with the model 
of Dovciak et al. (2004) to confront the issue of determining the spin of the black hole 
in MCG-6-30-15. Because of the extremely robust, broad, well-studied iron line in this 
system, MCG-6-30-15 is an excellent candidate for such a study. Given the complexity of 
the spectrum displayed by the source, it is important to perform this exercise in a step- 
by-step manner, clearly enumerating all of the assumptions at each stage, and employing 
physical models for the spectral complexity whenever possible. 

This guides the study presented in this Section. Due to the unprecedented signal-to- 
noise, we use the EPIC-pn data from the aforementioned long XMM-Newton observation of 
MCG-6-30-15. Data preparation and reduction followed Vaughan & Fabian (2004) exactly. 
In this Section, we present a step-by-step analysis of these data using models of increasing 
complexity and physical realism. Initially, to illustrate the potential power of broad iron lines 
for spin determination, we modeled the 2 — 10 keV EPIC-pn as a simple power-law continuum 
modified by a broad iron line (and absorption by the Galactic column of 4.1 x 10 20 cm~ 2 
toward this source); this is comparable to the study performed by Dovciak et al. (2004). 
It does, however, neglect the significant effects that continuum curvature due to ionized 
absorption could have on the inferred iron line parameters (and hence inferred black hole 
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spin). To assess these effects, we next model the effects of multiple warm absorbers and 
dust on first the 2 — lOkeV spectrum and then the full 0.6 — lOkeV spectrum. In our most 
sophisticated spectral model, we describe the 0.6 — lOkeV band including multiple absorption 
components and augmenting the simple broadened iron line with a relativistically smeared 
ionized X-ray reflection spectrum from Ross & Fabian (2005). 

The best-fit parameters and error bars for each progressive model fit are shown in 
Table 1. Each fit and the corresponding A% 2 with respect to changes in a in each case are 
shown in Figs. 6-10. The unfolded spectrum and the best-fit model components for both 
the single kerrdisk case and the full ionized reflection spectrum convolved with kyconv 
are shown in Figs. 11-12. For the warm absorber tables described in Table 1 and §4.2, all 
abundances are frozen at the solar value. For all model components, the redshift is set to 
z = 0.008, the optically determined value for MCG-6-30-15 (Reynolds et al. 1997). In all 
of the fitting below, the inner disk radius contributing to the iron line emission (or X-ray 
reflection spectrum in the case of Model 5) was not allowed to be smaller than the radius of 
marginal stability. 

4.1. Simple Power-Law Continuum and Iron Line Across the 2 — lOkeV 

Spectrum 

Initially, we perform an analysis of the 2 — lOkeV spectrum assuming that the under- 
lying continuum is a simple power-law (absorbed by the Galactic hydrogen column) and 
that the disk spectrum is just a single iron line (rather than a whole reflection spectrum). 
Fig. 4 shows the hard spectrum as fit by this photoabsorbed power-law (Model 1). To ac- 
curately model the continuum we have initially ignored the 4 — 7keV range when fitting 
this component. This prevents any contamination of the fit by the presence of an iron line 
reflection signature. Once the fit was complete, energies from 4 — 7keV were included again. 
As suspected, the data/model ratio shown in the lower panel demonstrates a significant 
residual feature above the continuum, which appears to have the form of a highly broadened 
iron line peaking at 6.4 keV. The presence of this large residual feature results in a poor 
fit, x 2 /dof = 4577/1106. Fig. 5a plots the data/model ratio again, this time for a model 
including a single, broad kerrdisk line with a rest-frame energy of 6.4 keV (Model 2). Two 
narrow redshifted gaussians representing a cold Fe-Ka line and an ionized line of iron are 
also included at 6.4 and 6.9 keV, respectively, as in Fabian et al. (2002). The 6.4 and 6.9 keV 
lines both have equivalent widths of ~ 14 eV. 

It should be noted here that Fabian et al. (2002) acknowledged that a 6.74 keV absorp- 
tion line of EW = —138 ± 35 eV fit the data as well as an emission line at 6.9 keV and 
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EW — 18 ± 6 eV. This absorption feature was also preferentially used in the MCG-6-30- 
15 work of Vaughan & Fabian (2004), and was consistent with the prediction of Sako et 
al. (2003) based on an RGS observation of this source in 2001. However, the 6.7 keV ab- 
sorption line detected in the deep high-resolution Chandra/HETGS spectrum of Young et 
al. (2005) is significantly weaker than that fitted by Fabian et al. (2002), with equivalent 
widths of EW = -18% -13 ± 9, and - 25 ± 9eV during the average, low and high flux 
states of the source, respectively. Thus, the high-resolution Chandra spectrum does not sup- 
port a spectral model for the EPIC spectrum in which the complexity in the 6.6 — 7.0 keV 
range includes a very strong absorption line. Two possible loopholes in this argument are 
(1) an order of magnitude temporal change in the helium-like iron column density between 
the XMM-Newton and Chandra grating observations and (2) extreme ~ 10 4 kms _1 velocity 
broadening of the absorption feature which would diminish its detectability in the high- 
resolution spectrum. More plausible is the notion of a weak (EW ~ 40 eV) and slightly 
broad emission line from hydrogen like iron (Fexxvi). Given that Chandra shows there 
to be a weak (EW ~ — 20 eV) narrow Fexxvi absorption line at 6.97keV, the EPIC-pn 
spectrum would be expected to show a net emission feature with EW ~ 20 eV. 

Most importantly, however, the details of whether this spectrum complexity is described 
by an ionized iron emission or absorption line has almost negligible effect on the broad iron 
line. In terms of the effect on the overall fit, Fabian et al. (2002) note that the relativistic 
line parameters differ insignificantly when one employs an absorption line rather than an 
emission line in the model. We have checked this result in our own analysis; specifically, we 
have replaced the 6.9 keV gaussian in Model 4 with a 6.74 keV gaussian in absorption (i.e., 
negative flux). The result was a Ax 2 /Adof = +26/0, indicating a marginal decrease in the 
overall goodness of fit. Visually, this fit was indistinguishable from the best fit from Model 4, 
and as per Fabian et al. , we also found minimal change in the relativistic line parameters. 
The inner emissivity index of the disk became marginally steeper and its inclination angle 
increased very slightly, but the changes were well within the statistical error bars. The 
most interesting point, however, is that the best fit equivalent width of this absorption line 
was only — 21.3eV; much less than the — 138 eV found by Fabian et al. . Such a modest 
equivalent width in comparison to what should be necessary suggests that this line is not 
robustly wanted in our fit to the data. 

The three iron lines (two narrow and one broad) that we do choose to include signifi- 
cantly improve the fit (x 2 /dof = 960/1096), and succeed in modeling out the residual feature 
shown in Fig. 4. At 90% confidence for one interesting parameter, our best fit for Model 2 
indicates a very rapidly rotating black hole (a = 0.970^^11) ■ To gauge the sensitivity of this 
fit to the spin of the black hole, Fig. 5b plots the change in the goodness of fit parameter 
Ax 2 as a function of black hole spin parameter. This clearly demonstrates that x 2 improves 



-20- 



dramatically as one approaches very rapid spins. The equivalent width of the broad iron line 
in our best-fitting model is ~ 729 eV, quite a bit higher than the ~ 550 eV cited by Fabian 
et al. (2002), the 450 eV found by Ballantyne et al. (2003), or the ~ 250 eV found by Dovciak 
et al. (2004). This unphysically high equivalent width is clearly due to the simplicity of this 
model — the effects of absorption and the X-ray reflection continuum, in particular, will 
introduce additional curvature into the continuum, thereby alleviating the need for such a 
strong line. 



4.2. Modeling the Warm Absorber 

In order to accurately assess the width and morphology of the iron line in an AGN, 
it is imperative that the soft portion of the spectrum be modeled correctly. One must 
be concerned about confusing a broad red-wing of a relativistic iron line with continuum 
curvature resulting from a putative "warm absorber" present within the AGN system. In 
fact, it has been argued that broad iron lines may be entirely unnecessary once one correctly 
accounts for the effects of a warm absorber (see Sako et al. 2003 for a discussion of MCG- 
6-30-15; see Turner et al. 2005 for a discussion of NGC 3516). As mentioned in §2, a deep 
Chandra/ 'HET 'GS observation of MCG-6-30-15 fails to find the K-shell absorption lines of 
the intermediate charge states of iron predicted from a model in which a warm absorber is 
mimicking the whole relativistic red-wing of the iron line. However, the question remains as 
to the effects that warm absorption has on fitting of the extreme red-wing of the iron line 
which drives black hole spin constraints. Very long data sets are needed in order to obtain 
spectra with enough signal to put both broad line and warm absorber models to the test 
and pursue the question of their overlap with statistical validity. The ~ 350 ks observation 
of MCG-6-30-15 taken by Fabian et al. (2002), which we use here, is ideal for such a study 
because of its large number of counts and its resolution of the broad iron feature this galaxy 
is thought to possess. 

We have used the xstar spectral synthesis package for photoionized gases (version 
2.1kn3, Kallman 2005) to construct a grid of warm absorber models as a function of the 
absorbed column density Ah and ionization parameter £. The ionization parameter is given 
by the usual definition, 



where L\ is the luminosity above the hydrogen Lyman limit, n e is the electron number 
density of the plasma and r is the distance from the (point) source of ionizing luminosity. 
We constructed 20 x 20 grids of models uniformly sampling the (log Ah, log£) plane in the 
range A H = 10 20 — > 10 24 cm~ 2 and £ = 1 — > 10 4 ergcm _1 s _1 . While these were made to be 
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multiplicative models (i.e., absorber models that can be applied to any emission spectrum), 
the ionization balance was solved assuming a power-law ionizing spectrum with a photon 
index of T = 2. This is a good approximation to MCG-6-30-15. 

Initially, we side-step the complexities of the soft (< 2keV) spectrum and apply this 
warm absorber model to the 2 — lOkeV only. In many ways, neglecting any constraints 
from the spectrum below 2keV maximizes the impact that it may have on the broad iron 
line; the sole "job" of the absorption component in this setting is to attempt to fit the 
curvature of the spectrum otherwise attributed to the broad iron line. We do notice a 
modest reduction in the goodness of fit parameter compared with Model 2 (the simple 
power-law and kerrdisk model) A% 2 /Adof = —26/2. The warm absorber in the best 
fitting Model 3 is of modest optical depth and rather weakly ionized: iV H = 4.22 x 10 22 cm -2 
and log £ = 0.84. Although the change in the goodness of fit is not dramatic, the additional 
continuum curvature introduced by the warm absorber leads to a reduction in the equivalent 
width of the iron line from EW = 729 eV down to the more physically reasonable value 
EW = 521 eV. As can be seen from Table 1, the parameters that determine the best- 
fitting shape of the iron line (the emissivity indices, inner radius, break radius, outer radius, 
inclination of the disk and black hole spin) are essentially unaffected by the inclusion of 
this warm absorbing component. As shown in Fig. 6, a rapidly rotating black hole is still 
preferred in this case (a = 0.997^0^35). However, in contrast to the case with Model 2, the 
inclusion of a warm absorber to the 2 — lOkeV spectrum allows models with low black hole 
spin to fit the data adequately (with low-spin cases producing a goodness of fit parameter 
which is only A% 2 ~ 40 worse that high-spin cases) — in these cases, the 3 — 4 keV curvature 
is being modeled as the effects of warm absorption rather than the extreme red-wing of the 
iron line. 

Of course, the soft X-ray band (< 2keV) is extremely important for constraining the 
properties of warm absorbers; the opacity of most well-studied warm absorbers is dominated 
by oxygen and iron edge/line absorption in this band. Hence, to be complete, we must 
extend our study of the effects of warm absorption on X-ray reflection features to the full 
0.6 - lOkeV band. 

When fitting the 0.6 — lOkeV band, we find that we cannot produce an acceptable fit 
with a model consisting of a simple power-law and kerrdisk line subjected to the effects of 
Galactic absorption and a one-zone warm absorber. This is not surprising: it is generally 
thought that warm absorbers must be physically more complex than a one-zone model can 
account for, i.e., they cannot be well described by a single value of the column density and 
ionization parameter. Physically, the warm absorber likely represents an wind emanating 
from the accretion disk and/or cold torus surrounding the central engine and may well contain 
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dust. A continuum of ionization parameters is likely to exist along the line of sight. For 
computational purposes, however, it is convenient to approximate this as a discrete number of 
zones, each of which is characterized by a column density and ionization parameter. Adding 
in a second warm absorber dramatically improves the fit and, indeed, the two absorbers do 
seem to represent distinct "zones" of material based on their column densities and ionization 
parameters (see Table 1 for details). 

Even taking both warm absorber models into account, however, there still appears to 
be significant remaining absorption in the spectrum below 2keV, as well as a strong soft 
excess below ~ 0.7 keV. A strong edge due to the L3-edge of neutral iron (presumably in 
dust grains embedded within the warm absorber) has already been noted in high-resolution 
Chandra and XMM-Newton grating spectra of MCG-6-30-15 (Lee et al. 2002; Turner et 
al. 2003). Incorporating this edge into our fit (employing spectral tables kindly provided to 
us by Julia Lee) makes a significant visual and statistical improvement, largely explaining the 
unmodeled absorption mentioned above. It is worth noting that the edge demands quite a 
high column density of iron: log(Np c ) = 17.54, which is approximately a factor of two higher 
than that found by Turner et al. (2003). To address the soft excess seen below 0.7 keV, 
we employ a simple blackbody model. Since the data are only sensitive to the tail of this 
component, the other parameters describing the spectrum are completely insensitive to the 
precise model used for this soft excess. Table 1 details the fit, which includes two warm 
absorption zones, the Fe-L 3 edge and the additional soft excess. 

Although we cannot statistically compare the Model 3 and Model 4 fits (since we have 
expanded the energy range of study between these models), Model 4 does appear to describe 
the full 0.6 — lOkeV spectrum very well (see Fig. 7). As before, the data strongly prefer a 
rapidly rotating black hole. The inclusion of the 0.6 — 2keV data apply extra constraints 
on the warm absorbers; the partial degeneracy found in Model 3 between the red-wing of 
the iron line and the curvature introduced by warm absorption is now removed. At 90% 
confidence, this model gives formal constraints on the black hole spin of a = 0.997 ± 0.001, 
and the broad iron line has an equivalent width of ~ 926 eV. This is significantly broader 
than we find for Models 2-3, or in any of the other analyses of this data set, and reflects 
the breadth of the red-wing of the iron line feature in this fit. The values for the best-fit 
parameters for Model 4 are shown in Table 1. 

Fig. 9 shows the unfolded spectrum for MCG-6-30-15 fit with Model 4 using a simple 
kerrdisk line. Each model component is colored and labeled separately to highlight its 
relative contribution to the fit. Note the relatively strong blackbody component that must 
be included at ~ 0.1 keV in order to accurately model the spectrum below ~ 2keV, as well 
as the redshifted gaussian emission lines at 6.4 and 6.9 keV that must be added to the broad 
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neutral iron line to fully capture the shape of the hard spectrum. 

Turner et al. (2003) have approached the question of absoprtion in MCG-6-30-15 by 
analyzing the RGS spectrum from the same XMM-Newton observation used here (from 
Fabian et al. , 2002). In fitting the 0.32 — 1.7keV range, these authors have identified six 
components of absorption: absorption by the cold Galactic column, four "zones" of warm 
absorbing plasma, and an absorption L 3 edge of neutral iron (from dust embedded in one 
or more of the warm absorbing zones). Considering this detail in structure identified by 
Turner et al. , it might appear that our spectral model (which only requires two warm 
absorber zones plus the neutral iron-L3 edge to describe the spectrum) is inconsistent with 
the picture painted by the RGS. Upon closer inspection, however, we see that this is not the 
case. Firstly, the lowest ionization warm absorber seen in the RGS (log£ ~ —4.42) cannot be 
distinguished from neutral absorption by our 0.6 — 10 keV EPIC-pn spectrum and hence is, 
in fact, accounted for through the neutral absorption column present in our model. Secondly, 
the two highest ionization warm absorbers identified by the RGS actually have rather similar 
ionization states (log£ ~ 1.6 — 1.7 for model-2 of Turner et al. , 2003) and are only separated 
into two zones through their kinematics; they are separated by ~ 2000 km s -1 in velocity 
space by the RGS. The EPIC-pn instrument, however, would not be able to resolve the 
velocity difference of these two zones. Accounting for these two facts, we would expect the 
Turner et al. four- zone RGS model to reduce to a two-zone model when applied to EPIC-pn 
data. 

The column densities and ionization parameters are lower in the Turner et al. fit than 
in ours, but unfortunately it is difficult to compare these values in a meaningful way due to 
calibration issues in the continuum response that presently exist with RGS data. This renders 
it nearly impossible to perform simultaneous RGS/EPIC-pn fits in XSPEC, which is why we 
have chosen not to address such a joint fit within the scope of this paper. Cross-calibration 
issues similarly affect our ability to compare EPIC-pn data with Chandra/ 'HET GS results, 
making it difficult to obtain a more precise, independent check on the fit to energies below 
~ 2 keV. 



4.3. Model Including a Full Reflection Spectrum 

The broadened iron emission line of Models 1-4 is, of course, just the tip of the ice- 
berg; the disk produces a whole spectrum of fluorescent and recombination lines, radiation- 
recombination continua and Compton backscattered continuum. To obtain truly reliable 
constraints, we must consider the full X-ray reflection spectrum. 
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In Model 5, we take the basic continuum/absorption components of Model 4 and aug- 
ment the simple iron line with the full X-ray reflection spectrum from an ionized disk surface 
(Ross & Fabian 2005). The Ross & Fabian (2005) models describe the reflected spectrum 
emitted by an optically-thick atmosphere (here, the surface of an accretion disk) of con- 
stant density that is illuminated by radiation with a power-law spectrum (here, photons 
that have been inverse Compton-scattered by relativistic electrons in the corona or base of 
a jet). We then convolve this reflection spectrum with the effects of relativistic smearing 
via kerrconv. Interestingly, as will be noted below, the soft X-ray emission associated with 
the photoionized disk surface naturally explains the soft excess without the need for an 
additional ad-hoc blackbody component. Hence, Model 5 does not include the blackbody 
component of Model 4. 

The best-fit parameters for this model are shown in Table 1, and in comparison with 
Models 1-4 it appears that Model 5 provides a statistical fit to the data that is not as 
good: (x 2 /dof = 1793/1374; A% 2 = +51 for one more degree of freedom compared with 
Model 4). Model 5 is, however, our most physical model in the sense that the whole reflection 
spectrum is treated as opposed to just the iron line, resulting in a natural explanation for 
the soft excess. That is to say that both the somewhat arbitrary blackbody component 
and the broad kerrdisk line of Model 4 are not required in this fit; the soft excess and 
broad iron feature are instead both fully described by the smeared radiative recombination 
line/continuum emission from the irradiated accretion disk. This change in the modeling 
of the soft excess increases the photon index of the continuum power-law component to 
T = 2.09, and also results in an increase in the inferred depth of the neutral iron-L 3 edge to 
an iron column density of log(A r F C ) = 17.68 (over 17.54 in Model 4). The column densities 
and ionization parameters of the two included warm absorbers also vary somewhat from 
Model 4, but the clear delineation between them remains evident. See Table 1 for details. 

The inclusion of an ionized X-ray reflection spectrum also has important implications for 
the derived spin parameter, as can be seen in Fig. 10. The fact that a significant component 
of the line broadening in this model is now due to Compton scattering reduces the inferred 
black hole spin to a = 0.9891q'qq2> slightly lower than the value determined in Model 4, but 
still consistent with a very rapidly spinning black hole. We do find, again, that a narrow 
Fe-Ka line at 6.4 keV is necessary in order to properly model the shape of the spectrum, as 
well as a 6.9 keV line of ionized iron as in Model 4. The equivalent widths of these narrow 
lines are 24.7eV and 27.3 eV, respectively. Including all of these components, we find that 
the total 0.6 - lOkeV luminosity of MCG-6-30-15 is L x = 9.34 x lO^ergs" 1 using WMAP 
cosmological parameters. 

Fig. 10 shows the plot of the relative contributions of the best-fit model components for 
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Model 5. The main features are an ionized disk reflection spectrum relativistically blurred 
with a kerrconv convolution model to represent the soft emission and ionized iron features, 
as well as two zgauss components to model the cold, neutral iron line at 6.4 keV and the 
6.9 keV line included in previous fits after Fabian et al. (2002). The absorption components 
in the soft spectrum are the same as those used for Model 4. Note that the presence of the 
ionized disk reflection negates the need for the blackbody component of Model 4 shown in 
Fig. 9. 

While not essential for the principal issue of this paper (i.e., determining black hole 
spin), it is instructive to estimate the "reflection fraction" of the ionized reflection spectrum, 
i? re fl. This parameter is defined to be proportional to the ratio of the normalization of the 
reflection spectrum to that of the intrinsic spectrum, and normalized such that -R rc fl = 1 
corresponds to a reflector that subtends half of the sky as seen from the X-ray source. 
Operationally, the ionized reflection model of Ross & Fabian (2005) is characterized by 
an absolute normalization and hence one cannot fit trivially for -R r efl- We estimate this 
parameter by extending Model 5 out to 100 keV and setting i? r efl to be the ratio of the 
normalization of the reflected and intrinsic spectra at the peak of the Compton reflection 
hump at 35keV. This technique yields _R rc fl ~ 1.25. Previous studies have also found the 
reflection parameter to be in this range for MCG-6-30-15, but it should be noted that in these 
cases i2 r efi = Vl/2ir has been a fitted parameter in the reflection model used, and as such has 
been considered an indication of the covering fraction of the reflecting material. For example, 
Fabian et al. (2002) began their spectral fitting by assuming i? rcfl = 1, but discovered that 
when this parameter was left free they achieved better statistical fits to the data. The best- 
fit value determined by these authors was _R re fl = 2.21q7, which is consistent with our own, 
within error bars. By contrast, Lee et al. (2000) used a ~ 400 ks simultaneous ASCAj RXTE 
observation of MCG-6-30-15 to identify four distinct spectral states for the source. These 
authors noted that as the flux increased, r 3 _ 10 steepened while r 10 _2o gradually flattened 
at approximately the same rate. Lee et al. modeled the reflection with a pexrav component 
in xspec, using a reflection inclination angle matching that of the accretion disk at % — 30° 
and possessing an exponential cutoff for the Comptonizing source power-law at 100 keV. The 
iron abundance was maintained at twice the solar value. These authors found that as the 
source flux increased, i? re fl did as well, beginning at 0.35^Q;}g and topping off at 1.371q;i4- 
The last two flux states are consistent with our result within error bars. 

A very high iron abundance in this disk is required by our Model 5 fit (Zp e > 9.8 Z Q 
at 90% confidence). Previous studies have also detected an overabundance of iron in the 
disk of MCG-6-30-15 (e.g. Ballantyne et al. 2003), though none have proposed this high 
an abundance. If we freeze the iron abundance to be 3 Z & and re- fit, as in Ballantyne et 
al. (2003), the goodness of fit parameter is increased by A% 2 = +463 for one less degree of 
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freedom, and there are obvious residual features created in the model fit by employing this 
tactic. Most of the absorption parameters remain relatively unchanged, but the disk param- 
eters are altered considerably: r min = 1.88 r g and r max = 12.26 r g , so the disk only radiates 
over a fairly thin ring near the radius of marginal stability. The emission profile consists of 
Oi\ = 5.60, «2 = 1-00, and = 3.22 r g , reinforcing this conclusion. The inclination angle of 
the disk is reduced to ~ 20° from ~ 30° in the best-fitting Model 5, and the spin of the black 
hole is considerably lowered to a = 0.31. Even though this is the best fit with Z Fe = 3 Z Q) it 
is clear from the form of the residuals around 6.4 keV that the reflection spectrum is being 
insufficiently broadened. Based on this fact and the substantial worsening of x 2 , it appears 
that our fit indeed prefers a higher abundance of iron than that found by Ballantyne et al. . 

Our confidence in the goodness of fit provided by Model 5 is strengthened by the consis- 
tency with which the continuum model we have employed extends to fit the BeppoSAX/PDS 
spectrum of MCG-6-30-15 at high energies. The BeppoSAX observation was taken simultane- 
ously with the XMM-Newton in 2001 and published by Fabian et al. (2002), acting as a valued 
"sanity check" to the derived EPIC-pn fit by providing spectral data from ~ 15 — 100 keV. 
We have performed a joint fit to the PDS and pn data within XSPEC, applying our Model 5 
to both sets of data. Even with the paucity of counts at high energies, the results of this fit 
clearly show that Model 5 provides an excellent description of both the pn and the PDS con- 
tinuum. The ability to place such additional constraints on our continuum and absorption 
parameters gives us greater confidence in the constraints we have correspondingly derived 
for the broad iron line parameters, including black hole spin. 



Table 1: Best-fit parameters for MCG-6-30-15 



Model Component 


Parameter 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


phabs 


AWcm -2 ) 


0.041 


0.041 


0.041 


0.041 


0.041 


UA 1 
Wii 1 
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i.yyo x iu 
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1.762 ± 0.085 


1.663 ± 0.033 
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/Vtjo f rm ~~ 2 ^ 

loP" ro 








3 489 X 1 23 
3.898 ± 0.17 
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flux (phcm~ 2 s _1 ) 
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1.907 ±0.007 


1.903 ±0.007 


1.966 ±0.014 


1.947 ±0.07 


2.090 ±0.007 




flux (phcm~ 2 s _1 ) 


1.387 ±0.009 x 10" 2 


1.363 ±0.01 x 10~ 2 


1.303 ±0.063 x 10" 3 


1.066 ±0.027 x 10" 4 


1.182 ±0.033 x 10~ 4 
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111 I ivC V ) 




a trja _i_ n n^7 

U.O(U HI Vl'KJOl 




(\ 4nn -i- n n^fi 
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kerrconv 


ai 




5.652 ±0.629 


5.461 ± 1.127 


6.555 ±0.428 


6.057 ±0.255 








2.863 ±0.151 


2.655 ±0.255 


2.439 ±0.265 


2.776 ±0.183 




n>r (r g ) 




5.747 ±0.816 


5.679 ± 1.622 


5.139 ±0.867 


5.555 ±0.662 




a 




0.9701°;™ 


n qq 7 +0.00l 
u - yy ' -0.035 


0.997 ±0.001 


n qsq+ 009 




i (°) 




20.405 ± 2.236 


23.800 ± 1.906 


29.700 ± 1.498 


29.699 ±0.709 




'"min (fg) 




1.738 


1.278 


1.278 


1.615 ±0.074 




^"max (?*g ) 




113.260 ± 1.692 


102.109 ± 1.178 


134.060 ± 0.638 


397.028 ±494.559 




flux ( ph cm S _1 ) 




2.980 ±0.359 x 10" 4 


1.781 ±0.332 x 10~ 5 


2.677 ±0.515 x 10" 6 




refl 


Fe / solar 

log £refl 
Trcfl 

flux (phcm~ 2 s _1 ) 










10.000 ±0.713 
2.027 ±0.001 
2.090 ± 0.006 
1.512 ±0.317 x 10" 7 


X 2 /dof 




4577/1106 (4.139) 


960/1096 (0.876) 


934/1094 (0.854) 


1742/1375 (1.267) 


1793/1374 (1.305) 



Models 1 — 3 are from 2 — lOkeV, Models 4 — 5 also include energies from 0.6 — 2.0 keV. Error 
bars quoted are all at the la level except for those on the spin parameter, which are all at 90% 
confidence. Errors for the column densities of the warm absorbers cannot be well constrained due 
to their magnitudes. For Model 5, the kerrconv component has no line energy or flux parameters, 
and the large error bars on the maximum disk radius indicates that this parameter is not robustly- 
constrained. 
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Fig. 4. — The phabs(po) fit to MCG-6-30-15. Notice the significant deviations of the data 
from this model, especially around 6.4keV. x 2 /dof = 4577/1106 (4.14). 
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Spin Parameter 



Fig. 5. — (a) A kerrdisk line near 6.4 keV has been added to the phabs(po) fit, as have two 
zgauss lines modeling out the narrow, cold iron line at 6.4 keV as well as the narrow, ionized 
iron line at 6.9 keV. x 2 /dof = 960/1096 (0.88). Note the flatness of the ratio plot shown 
here as compared with that shown above for the phabs(po) case, (b) The corresponding 
plot of the change in \ 2 vs - a f° r this model. At 90% confidence, a = 0.970io;oi5- 
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Fig. 6. — (a) A warm absorber table has been added to the fit to try and model the soft end 
of the spectrum more accurately. x 2 /dof = 934/1094 (0.85). (b) A% 2 vs. a plot. At 90% 
confidence, a = 0.997lg;^. 
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Spin Parameter 

Fig. 7. — (a) A second warm absorber has been added, as well as an iron edge at 0.707 keV 
and a blackbody component to model direct disk emission. The data now include energies 
from 0.6 — 2.0keV. The residual feature at ~ 1.8 keV is a calibration artifact (the Si K edge). 
X 2 /dof = 1742/1375 (1.27). (b) A X 2 vs. a plot. At 90% confidence, a = 0.997 ± 0.001. 
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Fig. 8. — (a) The bbody component has been replaced by a smeared ionized disk reflection 
spectrum. x 2 /dof = 1793/1374 (1.30). Again, the residual feature at ~ 1.8 keV is a calibra- 
tion artifact (the Si K edge). While this fit statistically seems less robust than that achieved 
with Model 4, the ionized reflection model is thought to be more physically accurate in its 
ability to account for both the soft excess in emission as well as the broad iron feature at 
6.4 keV. (b) A X 2 vs. a plot. At 90% confidence, a = 0.989ijj;££. 
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Fig. 9. — A vF v plot of the relative contributions of the model components for Model 4, shown 
above in Fig. 7a. The two redshifted gaussian components (zgauss) are shown in green and 
dark blue, the kerrdisk line is in light blue, the blackbody component (bbody) representing 
soft emission from the disk is in purple, and the power-law continuum (powerlaw) is in red. 
The soft components labeled "abs combo" in black represent the combination of absorption 
features present in the spectrum: Galactic photoabsorption, two warm absorber models and 
an iron absorption edge. The solid black line represents the combined model including all 
the dashed components. 
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Energy (keV) 

Fig. 10. — A vF v plot of the relative contributions of the model components for MCG-6-30- 
15, as above in Fig. 9. Now the fit is from Model 5, as shown in Fig. 8a. The color scheme 
is the same as for Fig. 11, but in this case the blackbody component has been replaced by 
an ionized reflection spectrum to model the soft emission as well as the broad neutral iron 
line. This component appears in light blue. 
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4.4. Ruling out a Schwarzchild Black Hole 

As mentioned above, Dovciak et al. (2004) and Beckwith & Done (2004) argue that 
broad lines cannot be used as black hole spin diagnostics due to the degeneracy that exists 
between the physical parameters that go into composing the line profile. We contend that 
this is not necessarily the case; broad lines can be used to constrain black hole spin if one 
takes into account the physical realism of the best-fit parameters. The degeneracy between 
parameters makes it difficult to calculate the precise angular momentum for a given black 
hole, but we can nonetheless statistically rule out certain regions of parameter space provided 
that the data used has the spectral resolution to enable accurate model fitting. 

In the case of MCG-6-30-15, the Fabian et al. (2002) data set is noteworthy for its length 
and unprecedented resolution of the broad iron feature. This makes it an ideal candidate 
for examining the parameter space of the new models in question. The width of the iron 
line implies that this feature is produced in the accretion disk immediately surrounding a 
rapidly spinning black hole, and indeed the best-fit kerrdisk parameters for the simplest 
model including an iron line (Model 2, from 2 — 10 keV) suggest that a = O^Olo^s with 
90% confidence (see Fig. 5a). The kerrdisk parameters are consistent with those of the 
kyrline best fit as well. Given that the fits excluding emission from within the radius of 
marginal stability imply a near-maximal spin for the black hole, here we pose the question 
"can we rule out the non-spinning case if we relax the restriction of no emission from within 
the radius of marginal stability?" . 

To answer this question we utilize the kyrline model, since kerrdisk does not yet 
include emission from within the ISCO. In previous fits to MCG-6-30-15 using low-spin 
black hole models, it has been found that the fit demands an inner emission radius well 
within the radius of marginal stability (Dovciak et al. 2004). Substituting two kyrline model 
components in for the kerrdisk component in the 2 — 10 keV fit, we freeze a = 0.0 and re-fit 
the data as in Model 2. Two kyrline components are used because the publicly released 
version of kyrline does not support a broken power-law emissivity index for the disk (as does 
kerrdisk), so we must divide the disk up into two effective regions: one extending outwards 
from r ms , and one interior to r ms representing the plunging region. Because the plunging 
region is physically distinct from the disk proper, any lack of continuity between either the 
emissivity indices or the fluxes of the two components is not considered problematic. 

The simplest kyrline fit for a Schwarzchild black hole is based on Model 2. Visually 
it does not differ perceptibly from the kerrdisk best fit for Model 2, and statistically it is 
a slightly better fit: A% 2 /Adof = —15/0 between the two fits. Most notable about the fit, 
however, is that when we force a = 0.0, the fit demands an inner radius that is deep within the 
plunging region (r min = 3.43 ± 0.19r g as compared with r ms = 6r g for a Schwarzchild black 
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hole) with an extremely high inner emissivity index a\ = 9.08 ± 1.36. Within the plunging 
region of a Schwarzschild black hole, we expect the radial component of the 4-velocity to be 



where, here, r is measured in units of r g = GM/c 2 (Reynolds & Begelman 1997). For 
r ps 3.43 we have u r ~ —0.22c, i.e., the material is already inflowing at mildly relativistic 
velocities. Hence, conservation of baryon number demands that this part of the accretion 
flow be extremely tenuous and, given that by assumption it is subjected to an intense X-ray 
irradiation, this material must be fully ionized (see discussion in Reynolds & Begelman 1997 
and Young, Ross & Fabian 1998). More quantitatively, the analysis of Reynolds & Begelman 
(1997; see Fig. 3 of this paper) show that the ionization parameter of the accreting matter at 
this radius will exceed log £ = 4 for any reasonable accretion efficiency; this is essentially fully 
ionized and will not imprint any obvious atomic signatures on the backscattered spectrum. 
It is therefore extremely unlikely that a disk with such a steep emissivity profile and an inner 
radius so deep within the plunging region is an accurate physical model of the real system. 

We have also performed a Schwarzchild fit to the data based on the more complex best- 
fitting model from 0.6 — lOkeV (Model 5). Recall that in this case we have introduced an 
ionized disk reflection component to the fit, which serves the dual purpose of modeling the 
soft excess and the broad Fe-Kct line at 6.4 keV. Whereas before in Model 5 we convolved 
our reflection spectrum with a kerrconv component, here we use a kyconv model instead 
because we are demanding that a = 0.0 in this case, and anticipate that a large fraction 
of the emission will come from within the ISCO, as was the case for our Schwarzchild fit 
based on Model 2 above. Following a similar approach outlined above, we use two kerrconv 
components to allow for a broken power-law emissivity index in the disk. Qualitatively 
similar results are obtained; the inner radius of the X-ray reflection is such a model is at 
r min — 3.02 ± 0.31 r g with an extremely steep inner emissivity index (a^ = 15.00 ± 3.12). 
As before with the simpler case of a Schwarzchild fit to Model 2, this implies that the vast 
majority of the emission of this model component originates well inside the ISCO, which 
is not physically realistic. The ionization parameter for the reflector has also risen for the 
Schwarzchild case from £ rcfl = 106.36 ± 0.94 to £ rcfl = 182.91 ± 13.58. The iron abundance 
has remained very high at Fe/solar = 10.0 ± 0.23, as was the case in Model 5. Even with 
these adjustments in parameters, however, the Schwarzchild fit visually fails to account for 
the entire breadth and shape of the 6.4 keV iron line. In this case x 2 /dof = 2855/1377 
(2.08), as compared with 1793/1374 (1.35) for the Model 5 fit where black hole spin is a free 
parameter. 

Based on these arguments for Model 2 and Model 5 (the best-fitting cases for 2 — 10 keV 
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and 0.6 — lOkeV, respectively) we can make a strong case that a non-rotating black hole 
cannot viably produce the broad iron feature in MCG-6-30-15. 

5. Discussion 

5.1. Summary of Results 

We have created a new model within xspec, called kerrdisk, which synthesizes fluo- 
rescent emission lines produced in the accretion disks surrounding black holes. Our model 
differs from the two most commonly used models (diskline and laor) in that it is fully 
relativistic and allows for the spin of the black hole to be a free parameter. We also enable 
the disk emissivity index to be modeled with greater precision as a broken power-law. While 
our model is comparable to the new models of Dovciak et al. (2004), Beckwith & Done (2004) 
and Cadez & Calvani (2005), we do not pre-calculate extensive tables of the photon transfer 
function in the manner of the first two sets of authors. Rather, we pre-tabulate a modest 
20 x 20 grid of slowly varying Cunningham transfer functions in (a, i) space over a given 
range of radii and redshifts within the disk, then linearly interpolate over these values to 
calculate the transfer function for a given point in (r, g*; a, i) space. The resulting line profile 
is indistinguishable from one produced by performing the relativistic calculations on the fly, 
in terms of accuracy, and is much faster. Making use of linear interpolation also negates 
the need for referring to gigabyte-sized tables (e.g. Dovciak et aland Beckwith & Done), 
making our code more portable. 

In fitting the hard spectrum of MCG-6-30-15 with the kerrdisk model, we have shown 
that the data prefer a fit with a spin parameter that tends towards the maximum value. 
A non-spinning black hole can produce a formally adequate fit (although still statistically 
worse than that achieved with a free spin parameter), but further requires a significant 
fraction of the X-ray reflection to originate unphysically deep within the plunging region. One 
might argue that for flows accreting at close to the Eddington rate, the radius of dynamical 
stability might be pushed to the marginally bound orbit, rather than the marginally stable 
orbit (Abramowicz et al. 1990; Chen et al. 1997). In principle this would mean that the 
optically thick part of the flow could come substantially closer to the event horizon, resulting 
in significantly more broad line contribution from this region. However, in the case of a 
Schwarzchild hole the marginally bound orbit is still at a radius of r m b = 4r g , which is well 
outside the minimum radius of emission we get for the Schwarzchild fit based on Model 5, 
where r min = 3.02 ± 0.31 r g . Also, it is unlikely that MCG-6-30-15 is accreting close to the 
Eddington rate. 
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In fact, for the rapidly-spinning black hole preferred by our spectral fitting, the issue of 
unmodeled emission from within the radius of marginal stability diminishes in importance. 
For such rapidly spinning black holes, the geometric area of the disk within this radius is 
small and the gravitational redshift of this region is extreme. Hence, we suspect that the 
inclusion of reasonable amounts of emission from within the radius of marginal stability will 
have negligible impact on our bet fitting spectral parameters, including the spin of the black 
hole. 

The best-fitting model for the full 0.6 — lOkeV data set appears to be Model 5. The 
hard spectrum is best represented by a power-law continuum and two gaussian features 
of iron at 6.4 (neutral) and 6.9 keV (ionized). The soft portion of the spectrum below 
~ 2 keV is well fit by relativistically blurred emission reflected from the surface of an ionized 
disk that is modified by a two-zone warm absorber, an iron absorption edge, and Galactic 
photoabsorption. This ionized emission also contributes to the breadth and morphology 
of the observed iron feature near 6.4 keV. As described above in §4.3, this model yields a 
best-fit spin parameter for the black hole in MCG-6-30-15 of a = 0.989+^^. Due to the 
aforementioned degeneracies in the broad iron line model parameters (§1) this value should 
not be interpreted as exact, but the fact that near maximal spin is approached in each of 
the Models 2-5 strongly indicates that the data prefer a rapidly spinning black hole. 



6. Conclusions and Future Work 

Both Dovciak et al. (2004) and Beckwith & Done (2004) make valid points about the 
difficulty of accurately calculating the spin parameter of a black hole based on broad line 
profiles from the disk alone. Given the number of parameters governing the kerrdisk line 
profile (nine, to be exact), some degeneracy between them is certainly to be expected, and it is 
possible to generate statistically indistinguishable line profiles using different combinations 
of parameters. This does not render the use of broad lines as spin diagnostics obsolete, 
however: provided that one uses care in determining whether parameter values are physically 
reasonable, it is possible to at least constrain the spin parameter to be within a certain range 
of values. In the example of MCG-6-30-15 above, we have clearly shown that the data rule 
out a non-spinning black hole if we also demand that the values for the disk emissivity indices 
are physically realistic. The plots of A% 2 vs. a clearly show the improvement in fit achieved 
when one frees the spin parameter. In this case the fit strongly tends toward maximal spin. 

The combination of such a high spin and such a high iron abundance in the disk may 
understandably give one pause when considering the precision of our best-fitting Model 5. 
Both a and Z Fe push the upper limits of their respective parameter spaces in order to try 
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to account for the extreme breadth and strength of the Fe-Kct feature in the spectrum. As 
we have mentioned, however, in this work we are seeking to establish a robust method for 
constraining black hole spin via X-ray spectroscopy and are not claiming that our best-fit 
parameter values should be interpreted as exact. The spin of the black hole may indeed 
reach a = 0.989, but there could also be some contribution to the line from highly redshifted 
radiation originating from within the ISCO that we have not yet accounted for, which could 
reduce the need for such a high spin value. What we can say with confidence is that fits 
to the data with a Schwarzchild black hole are not physically sound. Likewise, the disk 
may indeed have nearly ten times the solar value of iron, but underabundances in some 
of the lighter elements may also simulate an effectively iron-rich environment which could 
contribute to the observed width of the Fe-Ka line (Reynolds, Fabian & Inoue 1995). Suzaku 
observations of this source will be invaluable for untangling the parameters of the ionized 
reflection model. Much work remains to be done on the subjects of modeling accretion disks 
and isolating black hole spins, and improvements in the former will undoubtedly help us 
place more accurate constraints on the latter. 

If MCG-6-30-15 is indeed a rapidly spinning black hole, as seems likely given our spec- 
tral modeling, it is astrophysically interesting for several reasons. As mentioned in §1, rapidly 
spinning holes can, in principle, experience a magnetic torque by the fields threading the 
accretion disk at the radius of marginal stability. This torque can theoretically extract ro- 
tational energy from the hole itself, significantly enhancing the amount of dissipation in 
the inner accretion disk. The steepest dissipation profiles would be obtained if the mag- 
netic torque is applied completely at the radius of marginal stability (Agol & Krolik 2000). 
Therefore, only for rapidly spinning black holes would one expect to observe such a steep 
dissipation profile: in this scenario r ms is dragged inward very close to the event horizon, 
so the torque is strongest here. Such an effect would manifest itself via strongly redshifted 
reflection features in the spectrum, since the strong dissipation very close to the event hori- 
zon would mean that most of the emission would originate from this region. Based on its 
own best-fit emission profile, MCG-6-30-15 may in fact be giving us a glimpse of this phe- 
nomenon at work (Wilms et al. 2001). It should be noted, however, that interpreting the 
detected reflection features in this way demands that little of the observed emission originate 
from the plunging region within r ms . Taken in this context, the relatively steep emissivity 
index we found for the ionized disk in our best-fitting Model 5 (a = 0.989, a\ = 6.06) is not 
unexpected, and may be indicative of this type of magnetic torquing. 
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